Trapped Bose-Einstein condensates with Planck-scale induced deformation of the 

energy-momentum dispersion relation 
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We show that harmonically trapped Bose-Einstein condensates can be used to constrain Planck- 
scale physics. In particular we prove that a Planck-scale induced deformation of the Minkowski 
energy-momentum dispersion relation SE ~ ^\mcp/2M p produces a shift in the condensation tem- 
perature T c of about AT C /T° ~ 10 -6 £i for typical laboratory conditions. Such a shift allows to 
bound the deformation parameter up to |£i| < 10 4 . Moreover we show that it is possible to enlarge 
AT C /T° and improve the bound on £i lowering the frequency of the harmonic trap. Finally we 
compare the Planck-scale induced shift in T c with similar effects due to interboson interactions and 
finite size effects. 



A general feature of quantum gravitational theories, 
see Ref.[l[ for an introductory reading, is that they im- 
ply a deformation of the standard Minkowski free-particle 
energy-momentum dispersion relation at Planck scale, as 
in the case of loop quantum gravity and noncommu- 
tative geometries |7H10| . Such a feature opens a possibil- 
ity to test experimentally quantum gravitational effects 
111 ], The free particle deformed dispersion relation can 
be written in a general form as 



E(p) = E {p) + SE(p,m,M f 



(1) 



where Eq(p) = \/p 2 c 2 + m 2 c A is the Minkowski disper- 
sion relation, p the particle momentum, m its rest mass, 
c the speed of light and M p the Planck mass. The ex- 
plicit form of SE(p, m, Mp) depends on the details of the 
quantum gravity model used. For simplicity the relation 
([T]) is assumed to be universal, i.e., is the same for all the 
elementary particles, including composite particles such 
as nucleons or atoms when internal degrees of freedom 
are negligible. Since one should restore the Minkowski 
dispersion relation well above Planck scale, one should 
also require that SE(p,m, Mp) — > when M p — > oo. 
Moreover one would like to preserve the interpretation 
of m as the particle rest mass, and therefore one should 
impose the additional condition 5E{p = 0, m, Mp) = 0. 

In some case the deformation 5E can be tested in 
the ultra-relativistic regime p S> mc 2 by use of astro- 
physical data. In the specific case of the deformation 
SE(p,m, Mp) — rjip 2 /2M p such data could be sensitive 
to a deformation parameter \r]x\ ~ 1 as discussed by many 
authors [p34l6j . We mention that a preliminary analysis 
of the Fermi Space Telescope data [17H2l| is currently 
underway. 

Very recently 22, 23| it was argued the possibility to 
constrain the functional form of ([T]) in the nonrelativistic 
regime using ultra-precise cold-atom-recoil-frcqucncy ex- 



periments. Quite generally, in the nonrelativistic regime 
one can use the following low momentum (p ~ 0) asymp- 
totic expansion of SE(p, m, Mp) 

1 
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SE(p,m,M P ) 



2M V 



£,imcp + &p +£ 3 



(2) 



with the real deformation parameter £i associated to the 
leading term, £2 to the next leading term, and £3 to the 
next-to-next leading term. 

It might be objected that ((2|) can be ruled out for 
macroscopic objects when £1 ~ 1. From ([2]) one has 
p 2 /2m < SE for p < po = £1 cm 2 /Mp so that the defor- 
mation SE dominates over the Minkowski kinetic term 
of all momenta up to po- Standard-model particles with 
to < 10~ 16 Afp makes SE dominate in the extreme NR 



limit p 



< 



Po 



£1 10 



-16„ 



However, for macroscopic 



objects one can easily have m ~ Mp and po ~ £1 mc and 
therefore the deformation SE dominates over Minkowski 
kinetic term in the entire NR regime. But this would con- 
tradict the familiar dynamics of classical NR bodies and 
therefore rule out lp|). This is what is commonly named 
"soccer ball" problem, see 3, 2^. Note, however, that 
@ is merely the p ~ asymptotic expansion of the full 
deformation 5E(p,m, Mp), and is thus valid for all mo- 
menta up to some p\ , where p\ depends on the explicit 
functional form of SE. For example a deformation 



777c p 

SE(p,m,M Pl p x ) = tiTTTT- exp(-p/p x ) 



(3) 



< 



Px 



and SE ~ 



2M f 

behaves as SE — £ 1 mcp/2Mp for p 
for p^> p\. Therefore, to apply © to macroscopic bod- 
ies one should measure the momenta of extended objects 
with p < p\ and this is impossible for sufficiently small 
px below the lowest measurable momentum for extended 
bodies. Since one supposes that this is always the case 



2 



for pa, the relation ([2]) cannot be ruled out for classi- 
cal macroscopic bodies. We also emphasize how ([2]) is 
commonly accepted in the literature (22}-[23j]. 

In the preceding example p\ acts as a cutoff for the de- 
formation SE, above which Eq.® is no more valid. The 
value of the cutoff p\ depends on the specific quantum 
gravity model considered, but in principle it can be ar- 
bitrarily small. Therefore, in order to test the validity of 
Eq. (|5J| for a general class of quantum gravity models one 
must capture the effect of the deformation ([2]) at small 
momentums p < p\ for arbitrarily small p\. 

For that reason it is meaningful to look at the effect 
of the deformation ([2]) on the thermodynamic properties 
of Bose-Einstein consensates (BECs), which are, in the 
case of uniform free condensates, a collection of particles 
with p = momentum. Moreover, since in the case of 
harmonically trapped BECs the single particle ground 
energy is huj / 2 which corresponds to p — \j2timui where 
u) is the frequency of the harmonic trap, one also expects 
that Planck-scale effects are strongest for uniform (p = 0) 
than for trapped BECs. As we will show, this expectation 
is confirmed by calculations. 

In this Letter we calculate the shift in the critical tem- 
perature T c of a harmonically trapped BEC due to the 
Planck-scale induced deformation of the dispersion rela- 
tion @. The case of non-t rapp ed uniform BECs in a 
box has been studied in Ref.[26|. Neglecting interboson 
interactions and finite size effects in the BEC, it has been 
found that that the leading-order deformation in <j2j) pro- 
duces a shift 



AT C /T C ~ 0.1 & 



HM P n 1 / 3 



In [TV]. 



(4) 



where AT C = T c — T® with T c the critical temperature 
of the gas with the deformed dispersion relation, T° the 
condensation temperature in the undeformed Minkowski 
case, n the boson number density, m the boson mass, L 3 
the volume of the box and N = nL 3 the total number 
of particles. Such a shift is unexpectedly high if com- 
pared with the strength of the deformation in (J2]) which 
is of order SE/E ~ ^p/2cM P < £im/M p ~ 1CT 17 . 
For a condensate with N — 10 5 particles, a par- 

ticle number density n ~ 10 12 cm~ 3 and boson mass 
m ~ 150 x \Q- 27 kg one finds AT C /T C ° ~ 5.6 x 10~ 5 £i 
i26| but, since AT C /T C ° cx n -1 / 3 , such a shift can be even 
greater for dilute (smaller n) BECs. In fact, in the case 
of uniform BECs, one can reduce n and enlarge the tem- 
perature shift without however making n so small that it 
invalidates the thermodynamic limit. Alternatively one 
can increase the number of particles but, since the de- 
pendence on TV is logarithmic, one should consider huge 
N to considerably rise the temperature shift. Moreover, 
the effect of the next-to-leading-term deformation in ([2]) 
was found to be AT C /T® = ^Tn/Mp which is extremely 
small. 



The relevant result ((4]) encourages a further investiga- 
tion of the problem. In facts BECs are produced in the 
laboratory in laser-cooled, magnetically-trapped ultra- 
cold bosonic clouds (23, e.g. ^Rb 
HI, ifl_P, f 7 Rb MAHe & 41 
¥ Q 4 Yb dJand i<> [37|, 84 Sr[3 



2.: 
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34| J| 3 Cs IH, 
*Er jig. There- 
fore, to deal with laboratory measurements of AT C /T® 
one should generalize (jl]) to the case of trapped BECs. 

In the following we will show that in the case of har- 
monically trapped BECs one has a Planck-induced tem- 
perature shift AT C /T® ~ 10 -6 £i for the leading order 
deformation in ([2]). This allows to bound the deforma- 
tion parameter up to |£i| < 10 4 . Such a bound is four 
orders of magnitude above the best bound < 1 ob- 
tained with ultra-precise cold-atom-recoil-frequency ex- 
periments 



22|, |23[. However it is notable the possibil- 



ity to test Planck-scale effects with BECs. We will also 
discuss the possibility to improve such a bound consid- 
ering BECs with low values of the harmonic trap fre- 
quency uj. Moreover we will compare Planck-scale effects 
with the shift in the condensation temperature due to 
finite size effects and to interboson interactions. Finally 
we will show that the effect of the next to leading or- 
der term in Eq. ([2]) gives an extremely small temperature 
shift AT C /Tj? = ^m/Mp and therefore the deformation 
parameter £2 cannot be bounded significantly. 

Let us consider a system composed of N bosons 
trapped in an external spherically symmetric harmonic 
potential V ex t = muj 2 r 2 /2. We define <po(r) = 
(mijj /ttK) 3 / 4 exp[— V ex t(r)/ujh] that is the wave function 
of the ground state of the single particle quantum har- 
monic oscillator. The condensate distribution at zero 
temperature is that of N identical bosons in the ground 
state of a harmonic oscillator n c {r) = TV^M! 2 = 
N (y/nao)~ 3 exp[—mujr 2 /h} jj ^jj. The typical size of 
the support of ipo is ao = WK/muj and represents the 
size of the condensate at T — temperature. The 
typical size of the thermal cloud at finite t emp eratures 
k B T > htu is R T = a ^2Trk B T/huj > a [IJII3. We 
also remember that in the case of a harmonically trapped 
BEC the thermodynamic limit is given by N — > 00 and 
uj — » with Noj 3 finite. Finally we remember that 
the condensation temperature of the ideal harmonically 
trapped BEC with the usual Minkowski dispersion re- 
lation (£1 = £2 = £3 = 0) in the semiclssical limit 
k B T > hw is given by k B T° = Hoj (N/({3)) 1/3 > Huj. 

The semiclassical energy in phase space of nonrela- 
tivistic particles in presence of an external potential is 
E(p,r) = mc 2 + p 2 /2m + SE(p) + V ext {r), see H^fill- 
By means of this expression the number of bosons in 
thermal equilibrium Nth is given by 4oT43j] 



N th = 



d 3 xd 3 p 



(2nh) 3 exp[p(E(p,r) - /i)] - 1 ' (5) 
The expression §5§ allows to calculate the shift in the 
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condensation temperature due to the deformation ([2]) as 
follows. We parameterize the energy deformation ([2]) in 
the general form SE(p) = af(p) where a <C 1 is a dimen- 
sionless deformation parameter. With such a definition 
the condensation temperature T c (a) is a function of a. 
It is obtained by extracting T c (a) from Eq.([S]) with the 
substitutions Nth = N and /i = mc 2 , that is 



irh 3 N 



drdp r 2 p 2 



exp 



pV2m+q/(p) + V ext (r) 
k B T a (a) 



Since the lhs of ([6]) is independent of a one has d a N = 
and after some algebra one obtains 



d a T c (a) 
Tc(a) 



drdpr 2 p 2 f(p)g(p,r,a)p 2 / / drdpr 2 p 2 X 



[p 2 /2m + af(p) + V ext (r)] g{ Pl r, a) 



where 



g(p,r,a) 



exp 



p 2 /2m+q/(p) + y ext (r) 
t B T,(«) 



exp 



p 2 /2m+a f(p) + V s!ct (r) 



1 



(7) 



(8) 



We can use this expression to calculate the shift in T c 
Since a < 1, one has 



AT C _ T c (a) - T c (0) 



TP 



r c (o) 



d a r c (q) 
T c (a) 



(9) 



and the last term is evaluated by use of Eq.© remem- 
bering that T c (0) = T c °. 

We are first interested in the leading term of the defor- 
mation ((2J, that is SE — £imcp/2Mp which corresponds 
to a = £ > \m/2Mp <C 1 and f(p) — cp. In this case the 
resulting temperature shift is 



AT, 



/T c ° ~ 0.3 £i (m/M p ) mc 2 /hwN 1 / 3 . (10) 



and is finite in the thermodynamic limit, since it depends 
on uj and N through Nco 3 . For a 37-R& condensate in an 
harmonic spherical potential trap with typical values of 
N ~ 10 5 and frequency w — 10 Hz (iol - li^ . Eq.lQIj]) gives 



AT C /T C ° ~ 10~ 6 £i. This is about two order of magni- 
tude smaller than the non-trapped case ~ 10 _4 £i [26j j . 
but still extremely large if compared with the strength 
of deformation in Eq.® SE/E < 10~ 17 . This estima- 
tion of the Planck-induced AT C /T C ° should be compared 
with current experimental precision in T c measurements. 
In high-precision measurements of AT C /T® in 3 gK 44 1 
due to interboson interactions, the order of magnitude is 
AT C /T C ° ~5x 10~ 2 with a a ATc/T o < 10~ 2 . This allows 
to bound the deformation parameter up to |£i| < 10 4 , 
which is about four orders of magnitude above the best 



estimation J£i I < 1 obtained in cold-recoil- frequency ex- 
periments but is nevertheless notable. 

One can however look for physical situations in which 
the temperature shift is enlarged and the bound on £1 
improved. In Ref. 26] it was argued that in the case 
of uniform BECs the temperature shift ([4J is enlarged 
in dilute condensates since AT C /T® cx n^ 1 / 3 . Let us 
analyze the case of harmonically trapped condensates. 
The central density of an ideal trapped BEC is n c (0) = 
(6) N / (y/nao) 3 cx Nuj 3 ^ 2 and it goes to infinitum in the ther- 
modynamic limit. Moreover n c (0) is strongly suppressed 
in non-ideal BECs due to interactions (quantum deple- 
tion), so it is not useful to parameterize the tempera- 
ture shift. We can define the following number density 
n T = N/R 3 T = (Nlu 3 /c 3 )(mc 2 /2irk B T) 3 / 2 which is the 
ratio between the number of bosons and the volume oc- 
cupied by the thermal cloud and is well defined at finite 
temperatures in the thermodynamic limit. With such a 



-1/6 



therefore in the case 



definition one has AT C /T® cx 
of trapped condensates one can enlarge AT C /T° by reduc- 
ing riT which in turns implies to reduce Nlu 3 . Because 
of the —1/6 exponent the dependence of the temperature 
shift on tit is very smooth and it seems that one cannot 
significantly enlarge AT C /T° reducing tit- However one 
has ATc/T" cx lu~ 1 / 2 N~ 1 / 6 from which one see that one 
can enlarge AT C /T® lowering to. Also note that the de- 
pendence of the temperature shift with respect to N is 
very smooth and one cannot lower N without loosing the 
validity of the thermodynamic limit, therefore one can- 
not rise AT C /T® by lowering N. In conclusion in the case 
of harmonically trapped condensates one can rise Planck- 
scale effects on T c and therefore improve the bounds on 
£1 considering BECs with low tit, which corresponds to 
low frequencies lo and low condensation temperatures T c . 

We stress that an improvement of the precision in T c 
measurement would also allow to improve the bounds 
on £x- However, even in the ideal situation of an ex- 
tremely precise measurement of T c , in order to constrain 
Planck-scale effects with condensation temperature mea- 
surements one has to deal with the effect of interboson 
interactions and with finite size effects, which both affect 
T c . 

The shift in T c due to interboson interactions is 
AT C /T C ° ~ -3.426 a/X T where A T = h/ y/2irmk B T2 is 
the thermal wavelength at temperature and a is the 
s-wave scattering length 
as large as 4 x 10~ 2 for a 



4fj|-|42l|. Such a shift can be 
■ 10 2 a and A T ~ 10 4 a 
and it has been measured 



where anjs the Bohr radius 
in 3 gK [44[ where the data are well fitted by a second 
order polynomial AT C /T C ° ~ h{a/X T ) + 6 2 (a/A T ) 2 with 
61 = — 3.5±0.3 and 62 = 46±5, the second term being due 
to beyond-MF effects. Therefore in such a case the shift 
in T c due to interboson interaction is about four order 
of magnitude greater than the Planck-scale induced shift 
for £1 rv> 1. This problem may be solved in principle, just 
tuning the interaction coupling by Feshbach resonances 
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to very small values of a, therefore reducing the contri- 
bution of interboson interactions on AT C /T® below the 
Planck-scale induced shift ~ 10 -6 £i. However this may 
cause some practical difficulty. In fact, interboson inter- 
actions drive bosons to thermodynamic equilibrium and 
any experiment is limited by the fact that for extremely 
small interactions equilibrium is never reached. 

Finite size effects are due to the finiteness of N in com- 
parison with the thermodynamic limit N — > oo. The shift 
in the condensation temperature due to finite size effects 
is estimated as AT C /T C ° ~ -0.73 7V" 1 / 3 [H-g^ and typi- 
cally it is of about 1 -2 . Therefore, finite size effects are 
about four order of magnitude larger than the expected 
Planck-scale induced shift ~ 10~ 6 for £i ~ 1 and typi- 
cal N ~ 10 5 — I0 6 . We stress that finite size effects are 
technologically impossible to be tuned below Planck-scale 
induced effects at least for existing condensates, since 
this would require huge N > 10 30 corresponding to ex- 
tremely small shifts AT C /T C ° < 10~ 10 . Thus, in order to 
measure a deformation parameter |£i| one should predict 
the finite size contributions to AT C /T° for typical lab- 
oratory conditions N ~ I0 5 — I0 6 with an accuracy up 
to ~ I0 _6 |£i|. However, experimental strategies allows 
to eliminate finite size effects from AT C /T® measures as 



44j | . In fact for each measurement series at a given a 



and Xt, a reference measurement is taken with a small 
a/Xx ~ 0.005, same U) and very similar N, hence very 
similar Xt . Thus one can eliminate all a-independent sys- 
tematic errors that usually affect absolute measurements 
of T c including uncertainties in the absolute calibration 
of N and u as well as the shift due to finite size effects. 

We remark that, even if we have traced a possible ap- 
proach to the problem, dealing with interboson interac- 
tion and finite size effects in real experiments is a strong 
matter which deserves further investigation that will be 
presented elsewhere (iij . 

Finally we consider the next to leading order deforma- 
tion in Eq.® which is SE = ^2P 2 /2M p and corresponds 
to a = £,2m/Mp and f(p) = p 2 /2m. With such a de- 
formation Eq.® gives AT C /T C ° ~ ^m/Mp which is ex- 
tremely small of order ~ 10 -17 and does not allow to 
constrain significantly the next to leading order parame- 
ter £ 2 - 

In conclusion, the main goal of this Letter has been 
to show that trapped BECs can be used to constrain 
Planck-scale physics. In particular we have shown that 
the leading order dispersion relation deformation defined 
in Eq.@ produces a shift in the condensation temper- 
ature of about AT C /T® ~ 10 _6 £i for typical laboratory 
conditions and such a shift allows to bound the defor- 
mation parameter up to |£i| < I0 4 . Moreover we have 
discussed how it is possible to enlarge such a shift and 
improve the bound on £i lowering the frequency ui of the 
harmonic trap. Finally we have compared the Planck- 
scale induced shift with similar effects due to finite size 
and interboson interactions. 



During the edition of this Letter a similar analysis of 
the Planck-scale induced temperature shift in trapped 
BECS has been realized 0, However, in addi- 

tion to such analysis, here we give a stronger focus to 
the the possibility of planning specific experiments that 
might provide phcnomenological constraints on Planck- 
scale physics. 
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